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ABSTRACT 

We have assembled a spectroscopic sample of low-mass dwarfs observed as part of the Sloan Digital 
Sky Survey along one Galactic sightline, designed to investigate the observable properties of the thin 
and thick disks. This sample of ~ 7400 K and M stars also has measured ugriz photometry, proper 
motions, and radial velocities. We have computed UVW space motion distributions, and investigate 
their structure with respect to vertical distance from the Galactic Plane. We place constraints on the 
velocity dispersions of the thin and thick disks, using two-component Gaussian fits. We also compare 
these kinematic distributions to a leading Galactic model. Finally, we investigate other possible 
observable differences between the thin and thick disks, such as color, active fraction and metallicity. 

Subject headings: stars: low mass — stars: fundamental parameters — stars: M dwarfs — Galaxy: 
kinematics and dynamics — Galaxy: structure 



1. INTRODUCTION 

Modeling the Galaxy presents a challenging breadth 
of problems to both theorists and observers. Large N- 
Body simulations employing the constraints of gravity 
and A cold dark matter cosmology have sought to recre- 
ate the infant Galaxy, tracing the formation and col- 
lapse of baryons within the dark matter halo (Brook 
et al. 2004; Governatoet al. 2007). Observationally, these 
simulations are constrained by rotational velocities and 
luminosity profiles of extragalactic systems (Dalcanton 
et al. 1997). Closer to home, observers seek to recon- 
struct the history of the Milky Way's cannibalistic merg- 
ers through photometric identification of tidal debris, 
such as the Sagittarius dwarf (Ibata et al. 1994). Spec- 
troscopy is also employed to find surviving, co-moving 
stars of similar metallicity (e.g., Yanny et al. 2003). 

A particularly interesting problem to both the theorist 
and observer is the formation and nature of the thick 
disk (Gilmore & Reid 1983). This population has been 
explored extensively, mainly through star counts (e.g. 
Gilmore & Reid 1983; Reid & Majewski 1993; Buser 
et al. 1999; Norris 1999; Siegel et al. 2002). However, 
the thick disk scale height and local density normaliza- 
tion are still uncertain (Norris 1999), with values of /ithick 
ranging from ~ 700 - 1500 pc and local normalizations 
between 2% and 15%. Larger scale heights are usually 
coupled to lower normalization values (see Figure 1 of 
Siegel et al. 2002). The scale length of the thick disk is 
also uncertain, though typically scale lengths larger than 
the thin disk are inferred (Chen et al. 2001; Larsen & 
Humphreys 2003). The thick disk is thought to be an 
older, metal-poor population (Reid & Majewski 1993; 
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Chiba & Beers 2000). Metallicity differences, such as a 
element enhancement, have been explored (Bensby et al. 
2003), but other observables, such as photometric color 
and chromospheric activity trends, have yet to be stud- 
ied in detail. Thick disks observed in external galaxies 
(Yoachim & Dalcanton 2006) often appear to have struc- 
tural parameters and kinematics similar to the Milky 
Way. 

Modeling the observable properties of the thin and 
thick disks, namely star counts (Reid & Majewski 1993; 
Reid 1993) and kinematics (Mendez & van Altena 1996; 
Robin et al. 2003; Vallenari et al. 2006), has undergone a 
resurgence in recent years. The seminal work of Bahcall 
& Soneira (1980) laid the foundation for modeling star 
counts of the smooth Galactic components. Present-day 
models, such as the Besangon model (Robin et al. 2003), 
and the Padova model (Vallenari et al. 2006) also incor- 
porate kinematics, allowing for robust comparisons to 
observations. It is important to rigorously test their pre- 
dictions against actual observed spatial and kinematic 
distributions. 

Low-mass dwarfs are both ubiquitous and long-lived 
(Laughlin et al. 1997), and serve as excellent tracers of 
the Galactic potential (Wilson & Woolley 1970; Wielen 
1977; West et al. 2006). Modern surveys, such as the 
Sloan Digital Sky Survey (SDSS; York et al. 2000), are 
sensitive to K and early M dwarfs at distances of ~ 1 kpc 
above the Galactic Plane, probing the transition between 
the thin (Reid et al. 1997) and thick (Kerber et al. 2001) 
disks. At these distances, we expect about 20% of the ob- 
served stars to be thick disk members, assuming a local 
normalization of 2% (Reid & Majewski 1993) and scale 
heights of 300 pc and 1400 pc for the thin and thick disks, 
respectively. Star counts of low-mass dwarfs have been 
used to determine Galactic structural properties (Reid 
et al. 1997 and references therein). These studies sought 
to determine the vertical scale height and local normal- 
ization of the thin and thick disks (Siegel et al. 2002), as 
well as the underlying mass function of these populations 
(Martini & Osmer 1998; Phleps et al. 2000; Covey et al. 
2007a). The chemical evolution of the Galaxy has also 
been explored with red dwarfs (Reid et al. 1997), using 
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molecular band indices as a proxy for metallicity (Gizis 
1997). 

In addition to their utility as Galactic tracers, low-mass 
dwarfs have intrinsic properties, such as chromosphcric 
activity (West et al. 2004; Bochanski et al. 2005; Schmidt 
et al. 2007), that can be placed in a Galactic context. 
Additionally, metallicity may be probed using subdwarfs 
(Lcpinc et al. 2003), readily identified by their spectra, 
which show enhanced calcium hydride (CaH) absorption. 
Subdwarfs have been easily detected in large-scale sur- 
veys such as the SDSS (West et al. 2004). 

Kinematic studies of low-mass stars have a rich his- 
torical background. Samples are typically drawn from 
proper motion surveys, such as the New Luyten Two 
Tenths (NLTT) catalog (Luyten 1979) and the Lowell 
Proper Motion Survey (Giclas et al. 1971). Efforts have 
been made to identify nearby stars in these surveys, for 
example by Gliese & Jahreiss (1991). However, proper 
motion surveys possess inherent kinematic bias. The Mc- 
Cormick sample (Vyssotsky 1956), assembled from 875 
spectroscopically confirmed K and M dwarfs, has been 
frequently studied kinematically unbiased sample 

(Wielen 1977; Weis & Upgren 1995; Ratnatunga & Up- 
gren 1997), although it has been suggested that the sam- 
ple may be biased towards higher space motions (Rcid 
et al. 1995). The Palomar-Michigan State University sur- 
vey (PMSU; Reid et al. 1995; Hawley et al. 1996; Gizis 
et al. 2002; Reid et al. 2002) is among the largest prior 
spectroscopic surveys of low-mass stars, obtaining spec- 
tral types and radial velocities of ~ 1700 M dwarfs. The 
PMSU sample, which targeted objects from the Third 
Catalogue of Nearby Stars (Gliese & Jahreiss 1991), was 
used to construct a volume complete, kinematically un- 
biased catalog of ~ 500 stars, sampling distances to ~ 25 
pc. Other surveys, such as the 100 pc survey (Bochanski 
et al. 2005), have also used low-mass stars as kinematic 
probes. In Table 1, we summarize the sample sizes and 
approximate distance limits of previous major kinematic 
surveys of low-mass dwarfs, along with the mean velocity 
dispersions for each study. It is clear that our sample of 
several thousand stars out to distances of <~ 1 kpc re- 
sults in a study of Galactic kinematics using low-mass 
stars with unprecedented statistical significance. 

In this paper, we present our examination of the prop- 
erties of the thin and thick disks using a SDSS Low-Mass 
Spectroscopic Sample (SLoMaSS) of K and M dwarfs 
that is an order of magnitude larger than previous sam- 
ples. In §2, we describe the SDSS spectroscopic and pho- 
tometric observations that comprise SLoMaSS. The re- 
sulting distances and stellar velocities are presented in 
§3. In §4, we detail our efforts to separate the observed 
stars into two populations, search for kinematic, metal- 
licity and color gradients and compare our results to a 
contemporary Galactic model. Finally, §5 summarizes 
our findings. 

2. OBSERVATIONS 

2.1. SDSS Photometry 

The SDSS (York et al. 2000; Stoughton et al. 2002; 
Pier et al. 2003; Ivezic et al. 2004) is a large (-10,000 
sq. deg.), multi-color (ugriz; Fukugita et al. 1996; Gunn 
et al. 1998; Hogg et al. 2001; Smith et al. 2002; Tucker 
et al. 2006) photometric and spectroscopic survey cen- 
tered on the Northern Galactic Cap. The 2.5m telescope 



(Gunn et al. 2006), located at Apache Point Observatory 
scans the sky on great circles, as the camera (Gunn et al. 
1998) simultaneously images the sky in five bands to a 
faint limit of — 22.2 in r, with a typical uncertainty of 
~ 2% at r ~ 20 (Ivezic et al. 2003; Adclman-McCarthy 
et al. 2006). The last data release (DR5 6 ; Adelman- 
McCarthy et al. 2007) comprises 8000 sq. deg. of imag- 
ing, yielding photometry of — 217 million unique objects, 
including — 85 million stars. SDSS photometry has en- 
abled a myriad of studies on both Galactic structure (e.g. 
Yanny et al. 2000; Newberg et al. 2002; Juric et al. 2005; 
Belokurov et al. 2006) and low-mass stars (e.g. Haw- 
ley et al. 2002; Walkowicz et al. 2004; West et al. 2005; 
Davenport et al. 2006). 

2.2. SDSS Spectroscopy 

When sky conditions prohibit photometric observa- 
tions, the SDSS telescope is fitted with twin fiber-fed 
spectrographs. These instruments simultaneously obtain 
640 medium-resolution (R — 1,800), flux-calibrated, op- 
tical (3800-9200 A) spectra per 3° plate, permitting ra- 
dial velocity measurements for most stars with an un- 
certainty of — 10 km s _1 (Abazajian et al. 2004). A 
typical 45 minute observation yields a signal-to-noise ra- 
tio per pixel > 4 at g = 20.2 and i = 19.9 (Stoughton 
et al. 2002), with a broadband flux calibration uncer- 
tainty of - 4% (Abazajian et al. 2004). The DR5 sam- 
ple includes over 1 million spectra, with — 216,000 stellar 
spectra (Adelman-McCarthy et al. 2007). The majority 
of spectra in the SDSS database are drawn from 3 main 
samples, which target objects based on their photomet- 
ric colors and morphological properties. These samples 
are optimized to observe galaxies (Strauss et al. 2002), 
luminous red galaxies with z — 0.5 — 1.0 (Eiscnstein et al. 
2001), and high rcdshift quasars (Richards et al. 2002). 
However, low-mass stars have similar colors to some of 
these samples and therefore are observed serendipitously. 
SDSS spectroscopy of late-type dwarfs has been the focus 
of numerous studies (see Hawley et al. 2002; Raymond 
et al. 2003; West et al. 2004; Silvestri et al. 2006; West 
et al. 2006; Bochanski et al. 2007 and references therein). 

2.3. SDSS Low-Mass Spectroscopic Sample: SLoMaSS 

During Fall 2001, a call was placed to the SDSS col- 
laboration to design special spectroscopic plates that 
employed different targeting algorithms than the usual 
SDSS survey samples described in §2.2. We designed 
and proposed a series of observations to probe the local 
vertical structure of the Milky Way, obtaining spectra of 
low-mass dwarfs in the southern equatorial stripe (stripe 
82) of the SDSS photometric footprint during Fall 2002. 
This stripe at zero declination is repeatedly observed dur- 
ing the time of the year when the Northern Galactic Cap 
is not visible. These repeat scans sample an area of — 300 
sq. deg. and have been used to study Type la supernovae 
(Sako et al. 2005; Frieman et al. 2007), stellar variabil- 
ity (Sesar et al. 2007) and characterize the repeatabil- 
ity of SDSS photometry (Ivezic et al. 2007). SLoMaSS 
is comprised of stars with unsaturated griz photometry 
and extinction-corrected magnitude limits of 15 < i < 
18 and i — z > 0.2. A series of three spectroscopic tilings 

6 http://www.sdss.org/dr5/ 
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TABLE 1 

Previous Low-Mass Stellar Kinematic Survey Results 



Study 


Sample 
Size 


Distance 
Limit (pc) 


tu a 


c v 


cr w 




Comments 


Vyssotsky (1956); Reid et al. (1995) 


368 


~ 300 


39 


26 


23 McCormick stars in PMSU survey 


Wielen (1977) 


516 


~ 300 


39 


23 


20 McCormick stars 


Reid ct al. (1995) 


514 


25 


43 


31 


25 PMSU I volume-complete sample 


Ratnatunga & Upgrcn (1997) 


773 


~ 300 


30.6 


18.5 


7.4 McCormick stars 


Reid et al. (2002) 


436 


25 


37.9 


26.1 


20.5 PMSU IV volume-complete sample 




ZD 




1 s 

lo 


16 PMSU IV volume-complete sample, core 


Bochanski ct al. (2005) 


419 


100 


35 


21 


22 Non-active stars, core 


155 


100 


19 


15 


16 Active stars, core 


This study 


70 


100 


28.4 


21.3 


19.2 


z 


< 100 pc 






100 


25.7 


20.9 


14.1 


z 


< 100 pc, core 




4300 


500 


39.0 


30.0 


24.8 


z 


< 500 pc 






500 


34.1 


24.6 


20.9 


z 


< 500 pc, core 




6893 


1000 


44.1 


36.9 


27.9 


z 


< 1000 pc 






1000 


37.4 


27.2 


23.6 


z 


< 1000 pc, core 




7398 


1600 


46.6 


42.2 


29.3 Total Sample 






1600 


38.5 


28.6 


24.4 Total Sample, core 



The velocity dispersions in U ,V , and W arc measured in km s 



composed of 15 plates (numbers 1118-1132, centered on 
I ~ 105°, b ~ —62°) was observed, yielding a total of 
8880 stellar spectroscopic targets in SLoMaSS, each with 
ugriz photometry. 

3. ANALYSIS 

3.1. Spectral Types 

The spectroscopic data were analyzed with the HAM- 
MER suite of software (Covey et al. 2007b) which mea- 
sures spectral type, Ha emission line properties, and var- 
ious spectral band indices. This pipeline automatically 
assigns a spectral type to each star, then allows the user 
to confirm and adjust the spectral type manually, if nec- 
essary. The accuracy of the final spectral types is ± 
1 subtype. After examining the entire sample by eye, 
we selected spectroscopically confirmed K and M dwarfs. 
This cut decreased the sample from 8880 to 8696 stars. 
The final spectral type distribution of SLoMaSS, after 
the additional cuts explained below, together with the 
average i — z color for each spectral type bin, is shown 
in Figure 1. 

3.2. Distances - Photometric Parallax 

Distances to the SLoMaSS stars were determined us- 
ing the photometric parallax relations described in West 
et al. (2005) and Davenport ct al. (2006). The West et al. 
(2005) relation was employed for stars with i — z redder 
than 0.37 and the Davenport et al. (2006) parallax re- 
lation was applied to the bluer stars in SLoMaSS. We 
neglect reddening corrections, since the average extinc- 
tion computed for the total column along this line of sight 
from Schlegel et al. (1998) was small (< 0.05 in i), and 
most of these stars will lie in front of significant amounts 
of dust. Stars with i — z colors that did not fall within 
the appropriate boundaries of the photometric parallax 
relations (0.22 < i — z < 1.84) were removed from the 
sample, decreasing its size to 8280 stars. Additionally, we 
applied the photometric white dwarf-M dwarf pair cuts 
of Smolcic et al. (2004), removing 4 more stars from the 




K0 K1 K2 K3 K4 K5 K7 M0 M1 M2 M3 M4 M5 M6 M7 M8 
Spectral Type 




0.2 0.4 0.6 0.8 1.0 1.2 
/-z 



Fig. 1. — Upper panel: Distribution of spectral types within 
SLoMaSS. The mean i—z color at each spectral type is shown above 
each bin. Lower panels: i vs. i — z color-magnitude distribution and 
r — i vs. i — z color-color diagram of stars in SLoMaSS. The lowest 
contour (black) in the color-color diagram indicates a density of 
100 stars in a square 0.1 magnitudes wide. Each additional contour 
represents 100 additional stars. 



sample. The — z Hess diagram and r — i,i — z color- 
color diagram are shown in the lower panels of Figure 1. 
The distance from the Galactic Plane was computed as- 
suming the Sun's vertical position to be 15 pc above the 
Plane (Cohen 1995; Ng et al. 1997; Binney et al. 1997). 
The vertical distribution of stars in SLoMaSS is shown 
in Figure 2. Note that the decline in stellar density at 
large distances reflects the incompleteness of our sample, 
not the intrinsic stellar density, while the decline at small 
distances from the Plane is due to the saturation of SDSS 
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Fig. 2. — The absolute vertical distance distribution of the sam- 
ple binned every 100 pc. The solid thick histogram is the total dis- 
tribution, while the five remaining histograms represent the height 
distribution for five i — z color bins as described in the legend. Note 
that there are hundreds of stars per bin out to an absolute vertical 
distance of ~ 1000 pc. 

3.3. Distances - Spectroscopic Parallax 

In addition to the photometric parallax relations men- 
tioned above, we also computed the distances to the M 
dwarfs in SLoMaSS using the spectroscopic parallax re- 
lations of Hawley et al. (2002). These distances served as 
a check on the photometric parallax relations and were 
later used to divide the sample in a search for color dif- 
ferences between the thin and thick disks (see §4.3.2 for 
details). The spectroscopic parallax of K stars was not 
computed, as no reliable calibrated relations were avail- 
able. 

3.4. Radial Velocities, Proper Motions & UVW 
Velocities 

Radial velocities were computed by cross-correlating 
the M dwarf stellar spectra against the low-mass tem- 
plate spectra of Bochanski et al. (2007) with the IRAF 7 
task fxcor. This cross-correlation routine is based on 
the method described in Tonry & Davis (1979). The 
measured velocities have external errors of ~ 4 km s _1 . 
The K dwarf radial velocities were obtained directly from 
the SDSS/Princeton ID spectral pipeline 8 , with typical 
uncertainties of 10 km s _1 . 

Proper motions were determined from the 
SDSS+USNO-B proper motion catalog (Munn et al. 
2004). The catalog is ~ 90% complete over the mag- 
nitude limits of our sample, with random errors of ~ 
3.5 mas yr _1 . We removed an additional 878 stars with 
poorly measured proper motions 9 , reducing our final 
sample size to 7398 stars. 

7 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

8 http://spectro.princeton.edu 

9 Specifically, we required a match between catalogs (match > 
0), detections in at least 4 plates in USNO-B (nfit > 5) , no other 
objects within 7 arcscconds, which is the resolution of the POSS 
plates (dist22 > 7; sec Kilic ct al. 2006), and small errors in the 
proper motion determination (sigRA < 1000 and sigDcc < 1000); 
see (Munn ct al. 2004) for details. 



Combining distances, proper motions and radial veloc- 
ities, we computed the UVW space motions of each star 
in SLoMaSS using the method of Johnson & Sodcrblom 
(1987). The velocities are computed in a right-handed 
coordinate system, with positive U velocity directed to- 
ward the Galactic center and corrected for the solar mo- 
tion (10, 5, 7 km s _1 ; Dehnen & Binney 1998) with re- 
spect to the local standard of rest. 

4. RESULTS 

The SLoMaSS observations were used to investigate 
kinematics, colors, chromospheric activity and metallic- 
ity in the thin and thick disk populations. 

We used two methods in our analysis. First, we exam- 
ined the velocity distributions with no assumptions re- 
garding the parent population of a given star (see §4.1). 
These kinematic results were then used to test a leading 
Galactic model (§4.2). We also used the space motion of 
each star to assign it to the thin or thick disk popula- 
tion, and investigated the color, metallicity and activity 
differences between the two populations (§4.3). 

4.1. Kinematics 

The sample was binned in 100 pc increments of vertical 
distance from the Galactic Plane and probability plots of 
UVW velocities were constructed (Lutz & Upgren 1980; 
Reid et al. 1995, 2002). These diagrams plot the cumu- 
lative probability distribution in units of the standard 
deviation of the distribution. Hence, a single Gaussian 
distribution will appear as a straight line with a slope 
corresponding to the standard deviation of the sample 
and the y-intcrcept equal to the median of the distribu- 
tion. Non-Gaussian distributions will significantly devi- 
ate from a straight line. An example of velocity distribu- 
tions and their probability plots arc shown in Figure 3. 
An advantage to this method is its immunity to outliers 
and to binning effects from poorly populated histograms. 

For each distance bin in SLoMaSS, the probability 
plots are well fit with two lines: a low-dispersion, kine- 
matically colder "core" component (< | la\) and a high- 
dispersion, "wing" component (> | lcr | ) . Shown in Fig- 
ure 4 is an example of our analysis, illustrating linear 
fits to the core and wing distributions. At larger Galac- 
tic heights, the wing component traces the in situ thick 
disk population. 

The resulting dispersions are shown as a function of 
distance from the Galactic Plane in Figure 5. The low- 
dispersion core component is well-behaved, smoothly in- 
creasing with increasing height. The high-dispersion 
wing component is subject to larger scatter, but gener- 
ally increases with height from the Galactic Plane. Our 
results are summarized in Table 1. In order to facilitate 
comparison to the previous results included in the Ta- 
ble, we report the dispersions for several distance bins, 
as well as the entire sample. 

4.2. Galactic Dynamical Models 

4.2.1. Besangon Model: Introduction 

Our sample is well suited to rigorously test current 
Galactic models. We chose the Besangon model (Robin 
et al. 2003) as a fiducial, since it simulates both star 
counts and kinematics. This model is constructed on 
several assumptions and empirical constraints, with the 
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Fig. 3. — Schematics of velocity distributions (left panels) and 
their corresponding probability plots (right panels). The upper row 
displays two Gaussian distributions, one with a standard deviation 
of ~ 50 km s _1 (solid line) and another with a standard deviation 
of ~ 100 km s~ 1 (dashed line). A straight line in the probability 
plot signifies an underlying Gaussian distribution and the slope of 
the line is a measure of the standard deviation of the distribution. 
Note that the slope of the dashed line is twice that of the solid line. 
In the bottom row, a low dispersion component (dashed line) and 
high dispersion component (dotted line) are summed to produce 
the solid histogram. While the dual nature of the distribution is 
not readily apparent in the solid histogram, it is easily detected in 
the probability plot. 



Fig. 4. — Illustrative example of our velocity dispersion analysis. 
Shown arc probability plots for U, V and W for the 1200 < z < 1300 
pc bin (shaded lines). The "core" component (dashed line) is fit 
between ± 1 a, while the "wing" component (dot-dashed line) fits 
the outer edges of the distribution. Note the strong asymmetric 
drift component manifested as a change in slope at negative V 
velocities. 



goal of reproducing the stellar content of the Milky Way. 
Four populations comprise the model: the thin and thick 
disks, bulge and halo. For each population, a star- 
formation history, initial mass function, density law and 
age are imposed. Additionally, an age-metallicity rela- 
tion is employed for each population, with a Gaussian 
dispersion about the mean metallicity of each component 
(see Tables 1 and 3 of Robin et al. 2003). 

Kincmatically, the thin disk is composed of 7 groups of 
different ages (from 0.0 to 10 Gyr), each being isother- 
mal except for the youngest (0.0 - 0.15 Gyr). The thick 
disk is composed of an 11 Gyr old population, with a ve- 
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Fig. 5. — Velocity Dispersion as a function of vertical distance 
from the plane in 100 pc bins. The open squares represent the 
kincmatically colder, core component and the crosses represent the 
high-dispersion, wing component. Errors are derived from the lin- 
ear fits to the probability plots. Note the smooth trend in the 
low-dispersion component as a function of height from the Plane. 



locity ellipsoid based on the measurements of Ojha et al. 
(1996, 1999). An age- velocity dispersion relation from 
Gomez et al. (1997) is imposed on each component, and 
the model is then allowed to self-consistently evolve to 
the present-day. This self-consistency is achieved with 
the method described in Bienayme et al. (1987). Stellar 
populations are formed according to their appropriate 
density laws and introduced with an initial velocity dis- 
persion. The mass density of the stars is summed in a 
column of unit volume centered on the Sun along with 
dark matter halo and interstellar material contributions, 
and the potential is computed using the Poisson equa- 
tion. Stars are then evolved using the collisionless Boltz- 
mann equation in this new potential, and redistributed 
in the z direction, as complete orbital evolution is not 
included in this model. The process is iterated until the 
potential and scale heights of the disk populations con- 
verge at the 1% level. The local-mass density, which was 
empirically determined by Creze et al. (1998), also im- 
poses a constraint on the initial Galactic potential. We 
note that imposing an age- velocity dispersion relation on 
the model dominates the kinematics, and that there are 
significant uncertainties in these relations. For example, 
the age-velocity dispersions used by Rocha-Pinto et al. 
(2004) differ by factors of ~ 2 for the oldest stars, com- 
pared to the relations in Gomez et al. (1997). These 
discrepancies are primarily derived from the difficulty in 
determining ages of field stars, with systematic differ- 
ences imposed by various methods (i.e. isochronc fitting 
vs. chromospheric ages). 

Intrinsic properties such as age, mass, luminosity, 
metallicity, position and velocity are available for each 
star in the simulation. Observable properties, such as 
colors, proper motions and radial velocities are also re- 
ported. To determine colors, the model uses the Lejeune 
et al. (1998) database, which employs adjusted stellar 
synthetic spectra that attempt to match empirical color- 
temperature relations. 

4.2.2. Comparison to SDSS 
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We queried the Besangon webpage , generating a 
suite of synthetic datasets along the appropriate Galac- 
tic sightline of SLoMaSS (I ~ 105°, b 62°) with 

proper magnitude limits and error characteristics. We in- 
cluded both the thin and thick disks in the model inputs. 
The bulge and halo components were excluded since SLo- 
MaSS points away from the bulge, and with a maximum 
distance of ~ 2000 pc, we expect halo contamination to 
be minimal. A total of 25 models were generated from 
identical input, in order to minimize Poisson noise. 

As a consistency test of the Besangon model, we com- 
pared the model star counts to those obtained with SDSS 
survey photometry in the SLoMaSS field. We queried the 
SDSS Catalog Archive Server 11 for good point spread 
photometry 12 of stars along the appropriate sightline. 
This query resulted in 46,730 stellar targets. In the 
models, the mean star count was 46,991 with a stan- 
dard deviation of 137. This excellent agreement should 
be viewed cautiously, as only rough magnitude cuts were 
imposed on the Besangon models, and we did not at- 
tempt to model observational problems that would affect 
SDSS star counts (Vanden Berk et al. 2005), such as cos- 
mic rays or diffraction spikes near bright stars. However, 
the test inspires confidence that the Besangon model ad- 
equately represents the observed SDSS star counts. 

4.2.3. Comparison to SLoMaSS 

Each of the 25 models was sub-sampled to reproduce 
the color and distance distributions observed in SLo- 
MaSS. The two-dimensional color-distance density dis- 
tribution was computed for SLoMaSS, and color-distance 
pairs were randomly drawn according to this normal- 
ized density distribution. If the color-distance pair was 
found in the Besangon model, then it was kept. This 
Monte-Carlo sampling continued until there were 7398 
color-distance pairs, matching the number of stars in 
SLoMaSS. This sampling forced each model to simulate 
the properties of the SLoMaSS observations when com- 
pared to the complete SDSS photometry . Since SDSS 
r — i colors are not available for the Besangon model, 
the SLoMaSS colors were transformed to R — I using the 
relations of Davenport et al. (2006). Shown in Figure 6 
are the R — I and distance distributions for one instance 
of the model, along with those from SLoMaSS. Thus, 
each model is "observed" in a manner consistent with the 
spectroscopic observations of the SDSS field photometry. 

4.2.4. Kinematic Comparisons 

Using the sub-sampled model data, we compared the 
kinematic predictions of the Besangon models to the 
UVW velocities of SLoMaSS. Using the method de- 
scribed in §4.1, we constructed two sets of probability 
plots for each model: one composed of the thin and thick 
disk stars and one measured solely for the thin disk. The 
reason for this separation was twofold. The first was 
to isolate any systematics between the models thin and 

10 http://physique.obs-besancon.fr/modclc/ 

11 http://cas.sdss.org/dr5/en/ 

12 Specifically, we required the following flags: a detection 
in BINNED1 and no EDGE, NOPROFILE, PEAKCEN- 
TER, NOTCHECKED, PSF_FLUX_INTERP, SATURATED, 
BAD.COUNTS-ERROR, DEBLEND.NOPEAK, CHILD, 
BLENDED, INTERP.CENTER or COSMIC-RAY flags set 
(Stoughton et al. 2002). 



thick disk predictions. Additionally, this allowed for di- 
rect testing comparison of the thin disk predictions, since 
most of these stars lie on the inner "core" region of the 
probability plot. That is, if the velocity predictions are 
correct for the thin disk, the overall slope of the isolated 
thin disk models should roughly match the "core" slope 
measured from SLoMaSS. 

In Figure 7, an illustrative example of this analysis is 
shown. The main effect of adding the thick disk compo- 
nent is to increase the slope of both the core and wing 
components. Additionally, the wing component is en- 
hanced, as seen in the W velocity distribution, which is 
expected from addition of a high dispersion population. 
It is clear from comparing to Figure 4 that the combi- 
nation of the thin and thick disk model predictions are 
necessary to simulate the structure seen in the SLoMaSS 
data. 

Following the method explained in §4.1 we measured 
the slopes of the "core" and "wing" components of the 
each model as a function of distance from the Galactic 
Plane. The results of this analysis, compared to the SLo- 
MaSS results are shown in Figures 8 and 9. The SLo- 
MaSS results (open squares), which are shown in Fig- 
ure 5, are compared to the average Besangon prediction 
(crosses) for the thin (upper panels) and thick disk (lower 
panels) . While the model does well in predicting general 
trends, there are clearly some systematics. The model 
predictions for the aw thin disk velocity dispersions are 
systematically low, suggesting there may be flaws in the 
method used to compute these motions. Additionally, 
the model underestimates the thick disk oy dispersions 
at large Galactic heights. As described above (§4.2), the 
assumed age-velocity dispersions relations, which dom- 
inate the predicted kinematic structure, are uncertain 
and may contribute to this disagreement. This specula- 
tion is supported by the results summarized in Table 1, 
which demonstrate that most previous surveys (as well 
as our own) have measured velocity dispersions signifi- 
cantly higher than those predicted from the model. The 
mean velocity distributions for both SLoMaSS and the 
Besangon model are shown in Figure 9. Again, the model 
in general performs well, but there are evident systematic 
differences. The mean V velocities predicted for the thin 
disk by the model are systematically high, and the thick 
disk V velocities diverge at large Galactic heights. The 
first difference may be attributed to our chosen solar mo- 
tion values (10, 5, 7 km s _1 ; Dehnen & Binney 1998). A 
larger adopted value of Vq , such as the classic value of 12 
km s _1 (Dclhayc 1965), would move the mean velocities 
towards agreement. The discrepancy in the thick disk V 
velocities is probably due to small number statistics, as 
seen in Figure 2. 

4.3. Differences between the Thin and Thick Disks 

In order to examine observable differences between M 
dwarf members of the thin and thick disks, we kinemat- 
ically separated the sample using the method of Bensby 
et al. (2003). This technique selects outliers in the wings 
of the three-dimensional Gaussian velocity distribution, 
and computes the probability of these stars belonging 
to the thick disk. The thin and thick disk populations 
were characterized by the velocity dispersions in Bensby 
et al. (2003). In order to minimize systematics, the UVW 
velocities were re-computed using distances determined 
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Fig. 6. — Shown are the R — I and distance distributions of the 
stars in SLoMaSS (solid line) and a sampled model (dashed line). 
Note the agreement between the two datasets, indicating that we 
are sampling each instance of the Besan§on model in a manner 
consistent with the spectroscopic observations in SLoMaSS. 




Fig. 7. — Probability Plots of the U, V and W velocities for 
pure thin disk (shaded line) and thin + thick models (solid line) 
in the 1200 < z < 1300 pc bin. This height bin corresponds to 
the SLoMaSS data shown in Figure 4. Note that when the thick 
disk is added to the distribution, the overall slope increases, and 
the wing components are enhanced, as seen in the W probability 
plot. The combined thin and thick disk velocity distributions are 
also necessary to match the two-component structure seen in the 
SLoMaSS probability plots. 



with the spectroscopic parallax relations of Hawley et al. 
(2002), as described above in §3.3. Thus, UVW should 
not a priori vary systematically with color. This also 
limits the analysis to the 6577 M dwarfs in the sample. 

4.3.1. Metallicity 

The exact formation mechanism of the thick disk is 
uncertain (see Majewski 1993 and references therein), 
but there is evidence that the mean metallicity of the 
thick disk is lower than that of the thin disk (Rcid 
& Majewski 1993; Chiba & Beers 2000; Bensby et al. 
2003). Additionally, differences in the a-element dis- 
tributions have been shown to be distinct (Fuhrmann 
1998; Feltzing et al. 2003; Bensby et al. 2003). While 
direct metallicity determinations of M dwarfs are diffi- 
cult (Woolf & Wallerstcin 2006), proxies of metallicity 
have been employed (Gizis 1997; Lepine et al. 2003; Bur- 



gasser & Kirkpatrick 2006). These previous studies have 
used the low-resolution molecular band indices (CaHl, 
CaH2, CaH3, and Ti05) defined in Reid et al. (1995) to 
roughly discriminate between solar-metallicity, subdwarf 
([m/H] ~-1.2) and extreme-subdwarf ([m/H] ~ -2) pop- 
ulations. Adapting the methods of Lepine et al. (2003) 
and Burgasser & Kirkpatrick (2006) we computed the ra- 
tio (CaH2 + CaH3) / Ti05 for each star in the sample, 
which varies such that a larger value indicates a higher 
metallicity. The mean ratio for each spectral type for the 
thin and thick disk populations is shown in Figure 10. 
The two populations do not separate within the errors, 
suggesting that the observed metallicity distributions do 
not differ greatly. 

4.3.2. Color Gradients 

A lower mean metallicity could also manifest itself as 
a color shift at a given spectral type. Specifically, West 
et al. (2004) showed that at a given r — i or i — z color, 
low metallicity subdwarfs ([m/H] ~ —1-2) were redder 
in g — r. After SLoMaSS was kinematically separated, 
the mean u — g,g — r,r — i and i — z colors of the thin and 
thick disk stars were computed at each spectral subtype, 
as shown in Figure 11. There are no significant color 
differences between the thin and thick disks at a given 
spectral type. This further suggests that the SLoMaSS 
stars are not probing a large spread in metallicity. 

4.3.3. Chromo spheric Activity 

Finally, if the thick disk is an older system (as sug- 
gested by its hotter kinematics) , it should possess a lower 
fraction of active stars (West et al. 2004, 2006). The 
chromospheric activity timescale varies with mass, such 
that higher mass stars lose their activity sooner (after ~ 
1 Gyr) than their low-mass counterparts (~ 10 Gyr). To 
observationally quantify chromospheric activity, the Ha 
equivalent width (EW) is measured for each M dwarf. 
We employed the technique introduced in West et al. 
(2004) and described in Bochanski et al. (2007), select- 
ing chromospherically active and inactive stars at each 
spectral type. To be classified as active, a star must 
have an EW of 1 A, and pass additional signal-to-noise 
and error tests described in West et al. (2004). Figure 12 
shows the active fraction of stars as a function of spectral 
type for both disk populations. While most early- type 
M dwarfs (M0-M3) lose their activity quite rapidly, the 
later types (M5) possess smaller active fractions in the 
thick disk, and exhibit the expected behavior of older 
systems. However, these results are only suggestive, as 
populations older than ~ 4 Gyr require M dwarfs with 
types later than M5 in order to be effective probes of the 
age (West et al. 2007). 

5. CONCLUSIONS 

The Milky Way (along with other spiral galaxies) is 
evidently a composite of a few major smooth compo- 
nents (the thin and thick disks and the halo) and many 
smaller structures, such as the tidal debris streams from 
the Sagittarius dwarf galaxy. Using spectra, proper mo- 
tions and photometry along one Galactic sight-line, we 
studied the kinematic and structural distributions of the 
smooth thin and thick disks as a function of vertical dis- 
tance from the Plane. We fit two-component Gaussians 
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Fig. 8. — Thin Disk (upper panels) and Thick Disk (lower panels) velocity dispersions for SLoMaSS data (open squares) and Besancon 
model predictions (crosses). The model reproduces general trends, such as increased dispersion with distance from the Galactic Plane. The 
thin disk W velocity dispersions are systematically low, as are the thick disk V dispersions at large heights. This is probably due to the 
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to the UVW distributions as a function of height, plac- 
ing new constraints on the thin and thick disk velocity 
dispersions. 

This sample was then employed to test the predictions 
of current Galactic models. The Besangon model was 
chosen for comparison, since it is widely accepted as a 
standard and models the kinematics of stellar popula- 
tions. We generated a suite of 25 models, and each model 
was sampled to be consistent with the colors and distance 



0.4 



0.3 



o 
< 

§ 0.2 
o 



0.1 



0.0 



□ Thin Disk 
x Thick Disk 



M0 M1 M2 M3 M4 
Spectral Type 



M5 
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spectral type for the thin (open squares) and thick (crosses) disk 
populations. An older population would show a lower active frac- 
tion, as suggested by the thick disk stars in the M5 bin. 



distribution of SLoMaSS. The bulk kinematic properties 
of the data and model were compared. They agree to 
~ 10 km s _1 , placing a limit on the accuracy of model 
predictions. However, aw is poorly fit by the model, 
suggesting that further investigation into modeling the 
kinematics of the thin and thick disks is necessary. In 
particular, the age-velocity dispersion plays an impor- 
tant role in the kinematic predictions, and a more exact 
definition of this relation is required. 

SLoMaSS was divided kinematically, assigning mem- 
bership of each star to cither the thin or thick disk. We 
inspected the two populations as functions of spectral 
type, searching for observable differences in the metal- 
licity, colors and chromospheric activity. While there 
was little observed difference between the metallicity and 
color distributions, the activity fraction distribution sug- 
gests that the thick disk displays an activity level con- 
sistent with an older population, although this result is 
purely suggestive and needs to be re-examined with later 
spectral types. The lack of a strong observable signal 
may have several causes. Primarily, kinematic separa- 
tion of populations is not perfect, as stars with extreme 
kinematics in one population (i.e. the thin disk) can 
masquerade as members of the other population. This 
would dilute observable differences, and in our case, the 
more numerous thin disk population may be polluting 
the thick disk sample, despite our best efforts to mini- 
mize this effect. Secondly, the intrinsic properties (i.e. 
metallicity) of the stars that compose the thin and thick 
disks are drawn from overlapping distributions. Stars 
that are at the edges of these distributions would also 
blur the distinction among observable properties, cloud- 
ing the best efforts of the observer. Further investiga- 
tions using stellar tracers along the entire main sequence 
should alleviate some of these issues. 
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